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It is known that frequency feedback demodulators can show in- 
stability in their response to step changes (mistiming) in input fre- 
quency. This work reports on some mathematical analyses of this 
phenomenon as described by differential equations arising from simple 
IF and feedback filters in the demodulator. These equations are studied 
for local and global stability by geometric or phase-plane analysis, by 
means of Lyapunov functions, and by the topological Poincare-Ben- 
dixson methods. A typical result is for the case of no feedback filter 
and one-pole baseband analog of the IF filter, and states in physical 
terms that if the mistiming is not too big, specifically if 

| mistiming \ < (half -power IF bandwidth) (1 + feedback gain) 

then solutions which are bounded away from zero amplitude approach 
the natural equilibrium point. Examples are given in which a suf- 
ficiently large mistuning makes the equilibrium point unstable. 

I. INTRODUCTION 

The frequency feedback (or frequency compression) demodulator 
for FM signals was proposed by J. G. Chaffee 1 in 1937. After some 
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twenty-five years, Chaffee's idea was found to have a particularly 
fitting application in the satellite communications experiments Echo 2 
and Telstar, 3 in which there was a high premium on detecting a low- 
power wide-band FM signal in noise. Nevertheless, since its invention, 
little progress has been made in the mathematical analysis of this cir- 
cuit. Approximate methods of analysis and synthesis have been pro- 
posed, and some of them experimentally verified as useful. 4 - 5 However, 
except for unpublished works by S. 0. Rice and T. R. Williams, the 
nonlinear character of the circuit away from equilibrium positions has 
not been considered. 

It is the aim of this paper to formulate briefly one of the problems 
arising in the analysis of the FM with feedback (FMFB) receiver, 
namely that of stability of its response to step changes in input fre- 
quency. We shall write equations describing this response and present 
results about local and global stability of solutions for simple cases. 

II. CIRCUIT DESCRIPTION 

The FMFB receiver has been extensively discussed in recent publi- 
cations, 45 so only a brief description of it is included here. Roughly 
speaking, the receiver is a conventional FM demodulator, with a local 
oscillator whose frequency is controlled linearly by the output of the 
detector. The object of this control is to reduce the index of modulation 
at the output of the mixer, so as to be able to use a narrower IF filter 
than in a conventional FM receiver, and thus to eliminate some of the 
noise accompanying the input signal. The action of the circuit is to 
follow the slowly varying frequency of an FM wave while looking at 
it through a moving narrow frequency "window." 

The circuit is closely related to the phase-locked oscillator, but it is 
distinguished from that device by having amplitude effects absent in 
the latter. Mathematically this distinction takes the form that in 
FMFB there is an amplitude variable for every phase variable, while 
in phase-lock these variables do not appear. Their presence critically 
affects and complicates circuit analysis: thus the simplest FMFB 
equation is in two dimensions, while the simplest phase- lock equation 
is the pendulum equation, in one. The FMFB receiver resembles the 
phase-locked oscillator in that both devices work by phase-locking 
onto an FM wave that varies slowly over a limited range; if this range 
is exceeded locking fails and oscillation can set in. This phenomenon 
is well-known in phase-locked oscillators ; in FMFB receivers a similar 
behavior has been described by L. H. Enloe. 4 It is to this stability prob- 
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lem that we address ourselves, endeavoring an analytical study of the 
stability of simple differential equations describing the mistiming of 
the incoming signal away from the normal carrier frequency. 

A typical result we prove states that if the mistiming w r/ is not too 
big, specifically, in physical terms, if (for a one-pole baseband analog 
of the IF filter) 

|w rf | < (half-power IF bandwidth) (1 + feedback gain), 

then, for the simplest receiver, solutions which are bounded away from 
zero amplitude approach the equilibrium or critical point. This and 
similar results are proved by using the Poincare-Bendixson theory, or 
with the help of Lyapunov functions. 

III. EQUATION'S FOR RECEIVER WITH IDEAL DETECTOR 

We shall write equations for the FMFB receiver (see Fig. 1) under 
the assumption that it contains an ideal frequency detector. That is, 
we assume that if the signal leaving the IF filter is a(t) cos (ul + 0(0), 
then the detector produces the output 0(0- Let the mixer input be 

x c cos ait — x„ sin uit, 

and let the mixer multiply this input by 

2 cos M + Mt)) 

where /3 (in practice and here > 0) is the feedback gain, and <p is the 
feedback signal. It is assumed that the IF filter is tuned to the difference 
frequency w = ai, — w 2 , and can be represented by an impulse response 
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Fig. 1 — FMFB receiver, block diagram. 
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of the form 2/(0 cos cot, with /(•) a baseband response such that f(t) = 
for t < 0. The sum («, + w 2 ) components of the mixer are essentially 
removed by the IF filter, and will be ignored. The difference (to, - w 2 ) 
components at the output of the mixer are 

cos cot{x e cos (3<p + x, sin p<p\ - sin iat[x, cos fa - x c sin ftp}. 

The response of the IF filter to these components has the form 



cos 



cot [ /(/ - u)[x e (u) cos Mi) + x.(u) sin fo(u)\ du 

Jo 

- s'mcot [ j(t - u){x,(u) cos/fy(u) - x c (u)Binfo(u)\ 

Jo 



du 



+ terms around 2co 

+ terms representing initial conditions. 

We shall assume that the passband of /(•) is small compared to 2o>, 
so that the components around 2co may be ignored as well. 

To complete the loop equations we must indicate how the feedback 
<p(-) is determined from the output of the IF filter. We set, for t > 

Ue (t) = rc (t) -]_/"/(/- u){x c (u) cospcp(u) + x.(u) tin fo(u)\ du, 

Jo 

Um ({) = r .(0 +[/(/- m){*.(^) cosMtt) - x e (u)amfo(u)\ du, 

Jo 

where r,.( • ) , r 8 ( • ) represent the effects of initial conditions in the filter 
at t = 0. Exclusive of the carrier, the angle modulation of the IF 
output, 

u c cos ut — u, sin cot, 
is just 

6 - tan -1 — , 

Me 

corresponding to an instantaneous frequency, 

. u r ii, — u.u c 
9 = -~2 , 

a 

where a = (i^ + i<")' /2 - This is the output of the ideal detector. 

The feedback frequency <p(-) controlling the voltage-controlled 
oscillator is obtained by filtering (?(•)• Thus 
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<p(t) = r(t) + f k(t - u)6(u) du, t > 

where fc(-) is the impulse response of the feedback filter, and r(-) 
represents the effect of initial conditions at t = 0. 

IV. DIFFERENTIAL EQUATION FOR THE SIMPLEST CASE 

When the baseband responses /(■) and A-(-) correspond to filters 
with rational transfer functions, the integral equations for u,.(-) and 
««(•) can be turned into differential equations in a well-known way. 
In the simplest case, when there is no feedback filter and /(•! corre- 
sponds to a (one-pole no-zero) filter with transfer /i/(\ + s), we ob- 
tain the equations 

u c = —\u e + n[x e cos (J8 + x, sin 00] 
u, = —Xu, + n[x, cos /30 — x e sin /30] 
& u c (x. cos /30 - x e sin /30) - u,(x € cos /30 + x, sin /30) 

«r + M. 

The introduction of polar coordinates >/,. = a cos 6, u H = a sin simpli- 
fies these equations to 

= £ (a;, cos (0 + 1)0 - .r c sin (0 + 1)0) ... 

a (1; 

a = -Xa + /x(-)'r cos (/3 + 1)0 + .r, sin (0 + 1)0). 

We first consider the stability of equations ( 1 1 when the input to the 
demodulator consists of the carrier cos « n t alone, with no signal. In 
this case x c = 1 , x a = 0, and the equations are 

*--*■•«<*» + »• (2) 

d = -Xo + n cos (/3 + 1)0. 

We recall that the critical points of a differential equation x = v(x) are 
the points x in the phase-space at which v(x) = 0. Those of the system (2) 
are then the points in the a, plane at which simultaneously = a = 0, 
namely, 

a = x ' j§ + 1 ' n an mte 6 er - 

Because of the periodic dependence of the right-hand side of (2) on 0, 
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it is possible and convenient to define £ = [p + 1)9, to write (2) as 
f = -(0 + Djsmf 



(3) 



d = — Xa + m cos f, 



and to consider only principal values of f, and thus only the critical 
point (m/a, 0) in the plane specified by the polar coordinates (a, f). 

Theorem 1: The equations (3) are globally asymptotically stable for 
all positive X and n; all solutions tend to the critical point n/\, in an 
exponential manner; f is monotone, and a is either monotone or has one 
minimum. 

Prooj: We start with an heuristic direct analysis of the trajectories. 
Consider in Fig. 2 the circle C in the a, f plane defined by a = 0, that is, 
a = n/\ cos f. With x = a cos f , y = a sin f we shall examine the 
directions of the trajectories of (3) at points on C. The equation of C is 

(^ ,:■ 

Since C is the locus a = 0, it is apparent that on C each trajectory 
(has a tangent that) is perpendicular to the radius from the origin. 
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BY EQUATIONS (3) 
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Fig. 2 — Phase plane for no mistiming. 
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By symmetry about the x-axis, we can restrict attention to y ^ 0. The 
slope of C at x, y is 

while the slope of the line through x, y perpendicular to the radius 
from the origin is 

V 
X X 



Since for < x < n/\, 

V _ 

.x '"7V \J " "\ x 



-2*11*-? < 2 ^-r s 



we find 



? - * > - M - 

X /* 



or since y > and < a: < ju/X 



* = J- U - & I 

d:r 2// \X 



-5 



Thus every trajectory is entering C on a = except at the origin and 
at the critical point, f is decreasing in f > 0. If a trajectory ever crosses 
C it can never again recross it and must approach the critical point; 
in this case a has a single minimum. If a trajectory never crosses C, 
it must simply slip into the critical point, because then both a and f 
decrease and are bounded below. 

These preliminaries lead us to define the Lyapunov function 

V = 2 VX " a C0S V + 2^ Sin f) 

= - (distance from a, f to critical point) 2 

Evidently V ^ 0, and V = only at p/\, 0. The rate of change of V 
along trajectories of (3) is 

V = ad — a - cos f + a 7 f .sin f 

X A 
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0+1 



= -2X7 -^ sin 2 f < 

A 

except at the critical point, where V - t = 0. It follows from Theorem II, 
p. 37 of Ref. 6, that the system (3) is globally asymptotically stable: 
all solutions tend exponentially to the critical point with reciprocal time 
constant 2X. When X = p, 2\ has the physical interpretation 

2\ = 2 X (half-power IF bandwidth). 

V. MISTUNINC! IN THE SIMPLEST CASE 

Let us assume that in equation ( 1 ) we have 

jr s = sin uj, x c = cos u)jt, 

corresponding to the "mistimed" carrier input cosUi + t»d)t, or to the 
constant modulating signal o> d . The equation ( 1 ) assumes the form 

d = **mfat -(jS+ 1)0) 
a 

a = _xa + m cos fat - (/3 + l)B), 
or with t= mt - (fi + 1)0, 



(4) 



a = — Xa + n cos f . 
The critical point of this system is determined by the conditions 

It is important to note that because of the possibility of going to low 
amplitudes there always exist critical points, regardless of the value of 
o> d . This situation is in sharp contrast with the phase-locked oscillator. 
For a filterless phase-locked oscillator the equation corresponding to 
(4) would be 

f = co,, - fj. sin f , 

which has no critical point if w d > /x. Thus in phase-lock there is 
usually a critical frequency deviation above which locking is impossi- 
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ble for lack of critical points, and below which it may or may not occur. 
In the FMFB receiver, though, the critical points always exist but, as 
we shall see later, they are not always stable. 

We determine the stability of the critical point (5) by the standard 
method of linearization. The matrix 



M±» msf <4L+i2 sinf 



J' 


da f 




d . 
if" 


d . 
da J 





-Msin f 



-X 



of partial derivatives, evaluated at the critical point, is the matrix A 
appropriate for the linearized system. The determinant of (si- A) turns 
out to be 



s 2 + sX(/3 + 2) + X 2 (/3 + 1) + 



o>,/ 



+ 1 



with roots in the left half-plane. Hence the critical point is stable; in 
a neighborhood of it the trajectories approach it. 

Because of the symmetry of the equations, there is no loss of gen- 
erality in assuming, as we do henceforth, that o> d < 0. This convention 
is used in Figs. 3 and 4. 

Although we have not proved it, it is natural (and we conjecture) 
that a separatrix lies between solutions which pass around the origin 
in the upper half of the a, £ plane, and solutions which just miss the 
origin as they go past it in the lower half of the plane. Roughly speak- 
ing, the former pick up an extra 2ir of phase before settling. This 
separatrix may even be a fan 7 of solutions each of which goes into the 
origin, although in all likelihood it consists of a single trajectory. This 
conclusion is supported by a heuristic low-amplitude analysis of equa- 
tions (4) suggested by J. A. Morrison. He writes (4) as the single 
equation 

d£ = aw d - /iQ3 + 1) sin f 
da na cos f — Xa 2 

and for |cos f | = 1 drops terms of order a 2 , obtaining 



rff 



d 



ixa cos f j- = au d — n(0 + 1) sin f = pa -psin f, 



da 
whence by integration from a„ to a 

sin f — 



da 



Uda _ (On) ( • ,. co d a 

n(fi + 2)-\a) \ Sm f " M (/3 + 2) 
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a = /*/X 




Fig. 3 — Phase plane for mistuning. 
This formula suggests that if a point (a„ , £,,) , on the circle 

v(£ + 2) _ ; . 



a = 



w«j 



sin f 



and near the origin, is on a trajectory then a nearby point on the circle 
is on the same trajectory. In other words, a trajectory going through 
the origin does so like the circle above, which is tangent to but outside 
the circle f = with equation 

a = a£+il 8in f . 

To avoid difficulties we shall consider only trajectories which are 
bounded away from the origin. 
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We now address ourselves to the global stability of the mistimed 
equations (4). Since the system is two dimensional it is possible to try 
to use the topological Poincare-Bendixson methods. 7 This is most con- 
veniently done by passing to rectangular coordinates x = a cos £, 
y = a sin £ again, and calculating the divergence. We rewrite (4) as 



We find 



-X.r + 



= v(x, y) 



y 



-2^2 + M(/3 + 1) -,-T 

nfixy 



2 — WrfJ/ 



dx 

dx 

dy = 
dy 

div v = 



-X - 


• ** p 


+ 


ff 


-X - 


- „Rr X - 
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f 


*** (x 2 


+ 


ff 


-2X 


Ufa 

X + 


7 





-2X - & cos f. 




Fig. 4 — Details of the curve A. 
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Lemma 1: In the portion of the fourth quadrant comprised by an 
arbitrary neighborhood of the origin there is always a curve K, joining the 
positive x-axis to the negative y-axis, such that on K trajectories of (4) 
cross K in the outward direction, i.e., out of the part cut off by K near the 
origin, and into the part separated by K from the origin. (See Fig. 3.) 

Proof: Let K consist of the circle a = e from the /-axis down to the 
point where this circle crosses ri = 0, i.e., until cos f = t\/n. The 
Cartesian coordinates of this point are 



Xe 2 



e * /a 2n 2 
M M 

From here let K continue to the y-axis at slope 1, i.e., let it consist of 
that portion of the line 

€ . /~2 2T2 ^JL 

y = x V/i — e X — — 

which is between its intercept y - x a on the ?/-axis and the circle a = 0. 

Now on a = e inside a = we have a > 0, so on the circular part 
of K all trajectories are entering a > e, even at x , y . At x , y the 
trajectory is actually tangent to a = e, but for small enough e it is 
pointed in the direction of increasing x and so there too it must enter 
a > e. 

On the linear part of K we want to verify that 



dv _ -Mz + ^-^^T? 
dx \ • .. _l ,.A y 2 

-\X + jU - «4/ + MP ^2 _j_ ^2 



< 1 



if e is small enough. This is true because on the linear part of K we have 
x — > 0, y — * 0, and 

cos f = > a . -* 

v.i- + y 

all monotonely and uniformly, as € -> 0; near e = the denominator 
of dy/dx is close to /i(fi + 1) for (a;, y) 6 # , so on the linear part of K the 
trajectories are moving in the direction of increasing x at a slope 
dy/dx < 1; hence they are crossing K in the direction of increasing 
amplitude. This proves the lemma. 

Theorem 2: If \ <a d | < X(0 + 1), then every trajectory of (4) that is 
bounded away from the origin approaches the critical point a = p/\ cos f , 
f = tan -1 w d /X(/3 + 1). 
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Proof: It suffices to consider only &> d < 0. All trajectories outside 
a = n/\ have a < 0, so it is enough to consider those starting inside, 
because the others get there eventually. Consider a path starting in 
^ f ^ 3ir/2, f ^ 0, a ^ n/X, and bounded away from the origin. 
Either it stays in this region forever, or it reaches the fourth quadrant, 
or else it moves into f > 0. If it stays there forever then there is a 
closed region, free of critical points and excluding the origin, in which 
it stays. By a result of Poincare (Ref. 7, p. 232), this closed region also 
contains a closed path y, of period say t. Then because y is closed, 
if f (0) is on 7, 



f(r) - f (0) = f f(/) dt = 0. 



But this is impossible since f < throughout the region in question. 

In the third quadrant a trajectory can cross f = only once. The 
argument just given also shows that no path bounded away from the 
origin can stay in the region w ^ f ^ 37r/2, f > 0, a ^ /a/A. Thus all 
paths starting in the region a 5= y./\ and bounded away from the origin 
reach the fourth quadrant. 

The inequality in the hypothesis implies that the circle f = inter- 
sects the circle a = n/X. Away from the origin we have f < for x = 0, 
y > and £ > for x - 0, y < 0. On < x ^ p/X, i/ = the tra- 
jectories enter the fourth quadrant intersected with a ^ n/\. Since a is 
nonincreasing on a = ,u/A, it follows from Lemma 1 that there is a 
region R with these properties: 

(i) R is closed. 

(ii) R C jo ^ /i/X} H fourth quadrant, 
(iii) ({a ^ m/M ^ fourth quadrant — R) is in an arbitrarily small 

neighborhood of the origin, 
(iv) R is entered by the path under consideration, 
(v) R is invariant, i.e. maps into itself under the motion. 

Indeed R can he chosen to be a 2-cell (homeomorph of a disc). We 
note that the divergence is negative throughout R. It follows from the 
criterium of Bendixson (Ref. 7, p. 238) that R contains no limit cycles 
nor even an oval going to and from a critical point. Since R contains 
only one critical point, it can contain no path-polygon. Thus two of the 
three alternatives in Bendixson's theorem (Ref. 7. p. 230) are ruled 
out, and all paths starting in R go to the critical point. Since we can 
associate a region like R with any trajectory bounded away from the 
origin, the theorem is proved. 
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We remark that when A = /* the condition of the theorem can be 
rendered in physical terms as 

| mistiming | ^ (half-power IF bandwidth) (1 + feedback gain). 

We next show that a result similar to Theorem 2 can be obtained by 
a Lyapunov function argument. 

Theorem 3: If 

2X03 + 1)* 



Wrf < 



& 



then every trajectory of (4) that is bounded away from the origin 
approaches the critical point a = /x/A, £ = tan -1 w d /\(fi + 1 1. 

Proof: Consider the scalar function V defined by 

27 = (Xa - m cos r) 2 + 03 + l) -2 (aa. rf - m(/3 + 1) sin f) 2 
- a 2 + OS + l)" 2 (af) 2 . 
We find 

V = — X(Xa — m cos f) 2 



+ 



0«d 



0+1) 

X 



5 (Xa - n cos f)(aco (/ - m(/3 + 1) sin f) 



"FfT (aw " _m(/3 + 1)sinr) ' 2 - 



2(0 + 1) : 



— V" is a quadratic form in a and <zf with determinant 

-0a> d X 

203 + l) 2 P + 1 

which is positive whenever | co,, | < 2X(/3 + 1) /0. 

Consider now a trajectory bounded away from the origin. It is 
clearly bounded, so it has a positive limiting set r + which is invariant 
and to which it tends. There is a constant k such that the trajectory is 
entirely contained in a bounded subregion fi of [V < k}. Hence r 1 C £2 
and V = on T + . Since a is bounded away from on the trajectory, 
it follows that a = and f = on T\ Thus the trajectory tends to the 
critical point. (This is a variant of the argument for Theorem VI, p. 58 
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of Ref. 6.) Again, the condition of the theorem is that | co d | not be 
too big, viz., 

| co,, | < 2 (half-power IF bandwidth) X (0 + l) 8/ 7& 

where /3 is the feedback gain. 

VI. ONE-POLE FILTER IN THE FEEDBACK 

After the filtcrless case considered so far, the next simplest model 
for the FMFB receiver would have one-pole, no-zero niters both as the 
baseband equivalent /(•) of the IF response, and as the response 
k(-) in the feedback. This is the simplest case that has appreciable 
practical import: the IF filter corresponds closely to the one-mesh 
design described by Giger and Chaffee (loc. cit., p. 1119 and Fig. 5, 
p. 1120) ; the feedback filter and the gain /? are a rudimentary version 
of the dc-and-baseband amplifier sketched by these authors (loc. cit. 
pp. 1121-22.) 

In this case the differential equations for the system are 

u r = — \u e + n(x c cos P<p -f- x, sin @<p) 

u, = — X/*» + m(#» cos 0<p — x c sin 0<p) 

6 -2/ • . \ d , _! u. 

9 = a Cu c u, — u,u r ) = — tan — 

at u e 

x = -yx +86 
(f = x 

with a = (u 2 c + u]) as before, and 5/(7 + s) the transfer function of 
the feedback filter. Upon setting £ = + (l<p these simplify to 

I = fix + _ (#« cos £ — x c sin £) 

d = — Xff + fi(x e cos£ + x, sin£) 

x = —yx H (x, cos £ — z,, sin £). 

When the modulating signal is a constant a> d , then &,(£) = sin w d t, 
x c (t) = cos w d f, and with f (0 = w,^ — £(£) the equations become 

f = cu d — #c sin f 
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a = —\a + M COS f 

x = —yx + — sin f. 

Q 



(6) 



Let us note heuristically and physically that if 8 = y -> °o then the 
feedback bandwidth goes to oo and we obtain the equations (4) of the 
simplest case. 

We start with a study of the stability of the critical point £„ , a., , x 

defined by 



f = tan" 



Wd 



Ml + 



08 



a = - cos f o 



(7) 



■I'd — 



Swd 



7 + 05 

The matrix A = {dfi/dXj) of partial derivatives evaluated at the criti- 
cal point is 

r a 

X 2 tan f o 



-X 

M COS f o 

— Msinfn -X 



-/3 







gX 2 tan To 

OA . 7 

M COS f o 

The determinant of (sI-A) is 

(s + X)((s + X)(s + 7) + jSSX) + 5 2 X 2 /3tan 2 To + X 2 tan 2 f (« + 7) 

- s 3 + (2X + 7 )s 2 + (X7 + 05X + X 2 tan 2 f )s 
+ X 2 (05 + tan 2 t (p8 + 7)) 

= s 3 + CI2S 2 + aiS + do • 
A necessary and sufficient condition for stability is that 

eti , a 2 , a > and a 2 a! > a . 
The first three are clearly true, and the last amounts to 



(2X + 7X2X7 + X 2 + \08) + 



V7 + /35/ 



> X 2 /35 + 



2 2 
7^j 

7 + /3S 
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This is symmetric in ±a><j , and is true for |*> d | small enough. It becomes 
false for large |w (/ | if 2\ < J38. If A = /i. and y = 8, then these numbers 
are the half-power bandwidths of the IF and feedback filter respec- 
tively, and we may say in physical terms that if 

, „ , n IF bandwidth 2X 

= feedback gain < 2 X = — 

feedback bandwidth 5 

then a very large mistiming cannot affect the local stability of the sys- 
tem, but if (S exceeds twice the ratio of IF to baseband widths then 
sufficiently large mistiming will make the system unstable. This result 
was first observed in an unpublished work (although with some errors) 
of T. R. Williams. 

The global stability of the equations (6) for the case with a one- 
pole in the feedback is a far more difficult topic than the local. Natur- 
ally, as more complicated filters are assumed for the IF and the feed- 
back, the dimension of the problem goes up, and the kind of geometric 
analysis we are using here becomes virtually impossible. In particular, 
the Poincare-Bendixson theory used earlier is already unavailable in 
three dimensions, and also there seems to be no ready way to prove 
the boundedness of solutions. Nevertheless some information can be 
obtained from the construction of a Lyapunov function for the case of 
no mistiming; all attempts to extend the method to the case of mistun- 
ing have failed. 

Theorem 4-' If <o d = (no mistiming) then every trajectory of (6) 
that is bounded away from the line a = approaches the critical point 
given by (7). 

Proof : Consider the scalar function V defined by 



X ' A + 7 + # A 5, X + 7+S/3 

V is certainly positive along any trajectory satisfying the hypothesis. 
We find after a lot of elementary calculus that 

W A + 7 + 5/3 

Since the trajectory assumed in the theorem is bounded away from 
the line a = 0, it is easy to see from the equations (6) that it is bounded. 
Thus the positive limiting set of this trajectory is a nonempty, compact 
invariant set T + , to which it tends. There is a constant k such that the 
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trajectory is eventually in a bounded subregion ft of {V < k). Thus 
r + Q and V = on r + . Consider now the largest invariant subset M of 
(7 = 0)nsi. Clearly T + C M. Thus the trajectory tends to M. On M, 
V = 0; hence since the trajectory is bounded away from a = 0, we see 
also that d = and f = on M. Now the equations d = 0, f = 
define a spiral curve C on the cylinder \a = \i cos f by the formula 

\f/ = -gtan f, 

and M is an invariant subregion of this curve, bounded away from 
a = (which C is not). On C the vector field defined by the equations 
must either vanish, or else must point in the +i£ direction, or else point 
in the — \f/ direction; this is because there can be no motion in the a, 
f plane on d = 0, f = 0. If either of the second two alternatives holds 
at a point of C, that point cannot belong to M, because the trajectory 
through it would move off C and M Q C. Hence M consists of C-points 
at which the. field vanishes, i.e., M = {critical point}. (Cf. Theorem VI, 
p. 58 of Ref. 6). 

Try as we might (and we tried many Vs) we have not succeeded in 
proving a version of Theorem 4 in which there was mistiming. If the 
same V is used with a d ?* as was used in Theorem 4, then it is no longer 
true that V < 0; thus the results we feel are there still elude proof. 

VII. COMPENSATED ATTENUATOR IN FEEDBACK 

In private communication, L. H. Enloe and B. R. Davis have sug- 
gested that probably the most important practical FMFB receiver is 
one having a single-pole (as the baseband equivalent of the) IF filter, 
and a single-pole— zero-feedback filter, i.e., one with transfer 

*_±_£. (8) 

s -\- b 

In the time domain this acts like a delta-function plus an exponential. 
From formula (6) we see that the right-hand side of the equation for { 
can be thought of as the output of a filter whose response is a delta- 
function plus an exponential. This suggests that the analysis in Section 6 
can be made to cover the filter transfer (8) as well as the one-pole, 
because the differential equations are such as naturally to supply the 
constant in (8) if it is not present. 
Writing the transfer function as 

s + a = c 



s+ b s + 
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with c = a — b, the differential equations for the system heeome 

u e = — \u r + n(x e cos (3<p + x, sin 0<p) 
v, = — \v, + n(x t cos (3</> — x c sin jSp) 

« -2/ . . \ d , -1 M 4 

= a "(w r M, — w,?/ r ) = Titan — 
a/ u c 

v = 6 + y 

y — —by + c6, 
with a = (w« + wj) 1/a again. We note that 

6 = - (x. cos (8<p + 0) - x t sin (B<p + 0)); 
a 

now we set £ = f3<p + 8 and simplify the equation to 

£ = l8i/ + ^ + ^ (x. cos £ - .r e sin £) 
a = — Xa + mC'V cos £ + x, sin £) 
y = -ty -| ( if cos £ - x e sin £). 

It 

With the modulating signal a constant w d , we have as for equation 
(6), xjt) = sin o>,it, x,(t) = cos w d f, and we can set £(t) = u> A t — £(t) 
to obtain the equations 

03+ 1)m sin f 

f = a> rf — £?/ 

a = — Xa + m cos f (9) 

?/ = -by + ^sin f. 

These equations have the same form as (6) except that c replaces 5, 
b replaces 7, and there is an extra (0 + 1) coefficient in the sin term of f . 
The critical point is at 

a„ = - cos f n 

A 

f. = tan- "" 



X(/3 + 1)^1 + f 
2/o = y 03 + 1) tan $■„ . 
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The matrix of partial derivatives evaluated at the critical point is 



-«■ + » » + '?££ 



— «sin f 



-X 



[c03 + 1)x -* + »££* -»j 

The appropriate polynomial is 
s 3 + s 2 (X(/3 + 1) + b + X) + s(\b((3 + 2) 

+ X 2 (/3 + 1)(1 + tan 2 f ) + cX0(0 + 1)) 

+ &03 + 1)X 2 tan 2 f () + c\ 2 (3(0 + 1) + c(/3 + D/3X 2 tan 2 To 

= s 3 + a 2 s 2 + <M + flo • 

A necessary and sufficient condition for local stability is then that 
a„ , a* , and a 2 all > 0, and that o*ai > ««. . It is clear that the first 
three conditions are met whenever c > 0, i.e., a > b. The case a = b 
is degenerate and reduces in dimension to the filterless case of Section 
4. Also, a> is always positive. However, since 



tan f = 



Wd 



X(/3 + l)(l + f ) 
we can write ai as 

a , = X6 + X 2 (/3 + 1) + XO + 1)(& + c0) + 



|2 l2 

0) d 



(6 + c/3) 2 (/3 + 1) 



If 



6 0-1- 



1 



&+ 1 



> X + /3a 



then the sum of the first three terms is negative, and a x < for \m\ 
sufficiently small. Similarly 



a n = cX 2 /3(/3 + 1) + 



| co d I 2 (6 + eg) 
1 + f Y<0 + 1) : 
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which is negative for any oj^ if b + eft < 0, and is also negative for 
|o) d | sufficiently small if c < 0. The case c < is physically a bit 
strange, because it is equivalent to having positive feedback in one 
loop of the feedback path; thus it is not surprising that in this case 
there can be instability even for w rf = 0. 

In the case c > only the condition a- 2 a } > a» is of concern and this 
is 

(X/3 + 2X + b)\b(0 + 2) + X 2 (/3 + 1)(1 + tan 2 f„) + cX/3(/3 + 1) 

> (6 + /3c) (jS + 1)X tan 2 f + cX 2 /3(/3 + 1), 
which simplifies somewhat to 

X 2 6(/3 2 + 7/3 + 9) 

+ X 3 (/3 + 1)03 + 2) + c\m + 1)(6 + X^ + X) + X6 2 (/3 + 2) 

> _ 1^1" (x -> (/3 + l)f3c _ X 3 ( ^ + i) (/3 + 2)). 

(* + !)<! + f) 

Again, this is symmetric in ±u>,i , and is true for jo» (/ [ small enough. It 
becomes false for |w d | large if 

/3c > X(/3 + 2). 

We can think of the feedback path in this example as consisting of 
two parallel branches whose sum is added, one being an amplifier with 
gain, 1, the other being a single-pole filter S with dc gain 1 and (half- 
power) bandwidth 6, in series with an amplifier A of gain c/b. We 
can then say in physical terms, assuming A = /x, that if 

a / to i ™ v IF bandwidth . , . 

3 < (fi + 2) X 77,—^;^ X gam of A 

then even a very large mistiming cannot affect local system stability, 
but if not, then a sufficiently large mistiming can. The fact that 
(/? + 2) appears as a factor on the right shows how much the zero in 
the compensating attenuator helps prevent instability, in agreement 
with what has been observed in practice by L. H. Enloe and B. R. 
Davis (private communication). 

The global stability of the equations (9) may he studied by the same 
methods as were used in Section 6 for that of equations (6), but this 
topic is not pursued further here. 
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